UNCLASSIFIED 


AD.?iLy® 


/<V  D^K^duced 

t 

the 

ARMED  SERVICES  'KHNICAL  INFORMATION  AGENCY 


ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  pua^pose 
other  than  in  connection  with  a  definitely  i*eiated 
goverrment  procurement  operation,  the  U.  S. 

Government  thereby  incurs  no  responsibility^^ 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  ojiy  vay 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  other¬ 
wise  as  in  any  manner  licensing  the  holder  oj.  gjiy 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


^1' 

OD 

r-0 

00 


h 

c 


r 

L  , 

c  • 
c. 


h- 

'*=ur 

I 


M 


H 

I  ■ 


IttKi 


miH* 

bi 


H'', 


«ll 


Hill 


MEMORANDUM  REPORT  NO.  1405 
MAY  1962 


A  METHOD  FOR  DETERMINING  ESTIMATES  OF  THE 
DYNAMIC  UNBALANCE  AND  THE  THRUST  MISALIGNMENT  OF 
SPIN-STABILIZED  ROCKETS 

W.  J.  Sacco 


Department  of  the  Army  Project  No.  503-06-002 
Ordnance  Management  Structure  Code  No.  5010. 11.812 

BALLISTIC  RESEARCH  LABORATORIES 


ABERDEEN  PROVING  GROUND,  MARYLAND 


ASTIA  AVAILABILITY  NOTICE 

Qualified  requestors  may  obtain  copies  of  this  report  from  ASTIA. 

I 


The  findings  in  this  repoi’t  are  not  to  be  construed 
as  an  official  Department  of  the  Army  position. 


BALLISTIC  RESEARCH  LABORATORIES 
MEMORANDUM  REPORT  NO.  l405 
MAY  1962 


A  METHOD  FOR  DETERMINING  ESTIMATES  OF  THE  DYNAMIC  UNBALANCE  AND  THE 
THRUST  MISALIGNMENT  OF  SPIN-STABILIZED  ROCKETS 


W.  J.  Sacco 


Computing  Laboratory- 


Department  of  the  Array  Project  No.  505-06-002 
Ordnance  Management  Structure  Code  No.  5OIO. 11.812 


ABERDEEN  PROVING  GROUND, 


MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


MEMORAHDUM  REPORT  NO.  l405 


WJSacco/ic 

Aberdeen  Proving  Ground,  Md. 
May  1962 


A  METHOD  FOR  DETERMINING  ESTIMATES  OF  THE  DYNAMIC  UNBALANCE  AND  THE 
THRUST  MISALIGNMENT  OF  SPIN-STABILIZED  ROCKETS 

ABSTRACT 

Solutions  to  the  equation  of  motion  of  a  spinner  artillery  rocket 
about  its  center  of  gravity  undei-  the  influence  of  thrust  misalignment 
and  dynamic  unbalance  are  obtained.  A  method  for  determining  effective 
values  of  dynamic  unbalance  and  thrust  misalignment  is  given. 
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INl'RODUCTION 


In  the  mathematical  theory  of  the  motion  of  spin-stabilized  rockets 
during  the  burning  phase,  the  mathematical  model  is  a  system  of  ordinary  __ 
linear  differential  equations  with  variable  coefficients  and  various 
forcing  terms.  Solution  of  the  differential  equations  involves  contri¬ 
butions  due  to  unit  values  of  both  the  initial  parameters  and  coefficients 
of  the  forcing  functions.  We  shall  be  concerned  here  only  with  the  angular 
motion  about  the  center  of  gravity  and  with  forcing  terms  involving  only 
the  thrust  misalignment  and  the  dynamic  unbalance. 

The  motion  of  a  spinner  rocket  is  similar  to  that  of  a  spinning  top 
or  gyroscope.  In  addition  to  the  spinning  motion  there  exists  a  precession 
and  nutation  about  the  center  of  gravity.  The  thrust  misalignment,  and 
the  dynamic  unbalance,  are  the  principal  contributors  to  the  transverse 
angular  motion  of  the  rocket.  In  fact,  in  the  second  part  of  this  report 
we  ohall  postulate,  for  an  idealized  situation,  that  and  are  the 
sole  contributors  to  the  transverse  angular  motion;  thereby,  enabling  us 
to  give  a  method  for  determining  the  approximate  average  values  of  L  and 

The  coordinate  system  which  we  shall  use  to  describe  the  motion  of 
the  spinning  rocket  is  shown  in  Figure  I.  OXYZ  is  a  right-handed  system 
of  axes  with  the  origin  0  at  the  center  of  gravity  of  the  projectile. 

The  reference  axis  OX  is  taken  in  the  initial  direction  of  the  rocket 
axis,  OY  upward  in  the  vertical  plane,  throu£^  OX  and  OZ  to  the  right  as 
viewed  from  the  rear  of  the  projectile. 

We  can  imagine  a  sphere  with  the  center  at  the  center  of  gravity  and 
whose  radius  is  equal  to  the  distance  from  the  center  of  gravity  to  the 
nose  of  the  projectile,  and  then  picture  a  plane  perpendicular  to  OX  and 
tangent  to  the  sphere.  The  nose  of  the  projectile  can  be  thou^t  of  as 
tracing  out  a  curve  in  this  plane,  and  we  let  the  quantities  ^  and  ^  be 
the  rectangular  coordinates  of  a  point  on  the  curve. 

We  can  treat  the  plane  as  a  complex  plane,  with  the  real  axis  parallel 
to  OY,  and  the  imaginary  axis  parallel  to  OZ,  and  the  angles  measured  clock¬ 
wise  from  the  vertical  as  viewed  from  the  rear  of  the  rocket.  In  this 
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representation  the  complex  number  ^  defined  by  /  =  +  i^^  specifies 

the  orientation  of  the  instantaneous  axis  with  respect  to  the  initial  axis 
of  the  rocket.  It  is  clear  that  |  ^  j  is  the  sine  of  the  angle  between 
the  two  directions.  We  shall  assume  that  |  ^  |  is  less  than  0.1  radians 
for  all  time  values  during  burning,  which  enables  us  to  use  the  approxi¬ 
mations  of  replacing  the  sine  of  the  angle  by  the  angle  expressed  in 
radians  and  the  cosine  by  1.  Experimentation  has  shown  that  these  approxi¬ 
mations  are  valid. 

During  the  burning  period  of  a  spinning  rocket  the  gas  flow  throu^ 
the  canted  nozzles  produces  the  forward  thrust  on  the  rocket  and  the 
torque  about  the  longitudinal  axis  which  produces  the  axial  spin. 

Generally  the  line  of  thrust  of  the  rocket  will  not  pass  through 
the  center  of  gravity  of  the  rocket  but  will  miss  it  by  a  distance  L^. 

Ihis  produces  a  torque  about  the  center  of  gravity  resulting  in  a  cross 
spin.  The  distance  is  the  so-called  the  thrust  misalignment. 

The  axis  corresponding  to  the  least  polar  moment  of  inertia  is  called 
the  principal  longitudinal  axis  of  inertia.  Ideally  it  should  coincide 
with  the  geometrical  axis  of  symmetry  of  the  rocket.  When  the  principal 
longitudinal  axis  of  inertia  of  the  spinning  rocket  makes  a  complex 
angle,  with  the  geometrical  axis,  the  resulting  effect  is  characterized 
by  a  moment  about  an  axis  perpendicular  to  the  rocket  axis.  The  angle 
is  called  the  dynamic  unbalance  angle. 

For  spin-stabilized  rockets  the  major  lateral  dispersion-producing 
factors  are  the  (l)  dynamic  unbalance,  (2)  linear  thrust  misalignment,  (5) 
initial  yaw,  (k)  the  Initial  transverse  angular  velocity,  and  (5)  variable 
winds.  For  the  Idealized  case  considered  here,  the  last  three  factors  are 
eliminated. 

The  equation  of  motion  considered  here  describes  the  angular  rotation 
of  the  rocket  about  an  instantaneous  transverse  axis  throu^  the  center 
of  gravity,  when  the  center  of  gravity  is  fixed  in  space.  The  equation 
which  is  derived  in  references  {[  H  ^  and  ^CH  ^  ,  has  the  form 

0  -  2iqcrt^}  =  f(t)  (1) 
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where 


i= 

0  =  complex  orientation  angle  , 

2q  =  ratio  of  eixial  and  transverse  moments  of  inertia, 
a  =  angxilar  acceleration  of  the  rocket  about  its  longitudinal  axis, 
t  =  time  measured  from  ignition, 

f(t)  =  resultant  of  cross-torques  due  to  linear  thrust  misalignment 
and  dynamic  unbalance,  i.e. 


and 


-TL  e 
c 


-(l-2q)0^  ia-(at)^  j 


2  ]  iat' 
e 


=  transverse  moment  of  inertia. 


(due  to  thrust  misalignment) 
(dynamic  unbalance). 


T  =  thrust  of  rocket  • 


In  these  representations,  and  are  complex  parameters  which 
characterize  the  linear  thrust  misalignment  and  the  dynamic  unbalance 
respectively. 

Equation  (l)  is  a  linear  differential  equation  with  variable  co¬ 
efficients.  For  our  purposes  here,  q,a,B,T,  |  | |  |  assumed 

to  be  constants. 
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SOLUTIONS  TO  THE  EQUATION  OF  MOTION 


We  perform  two  integrations  on  equation  (l),  assuming  initially 
that  =  0,  and  obtain  the  expression 


|(t)  = 

t  t 

o  o 


e  dudx 


Case  a)  The  thrust  misalignment. 

We  let  ^^^(t)  equal  the  f(t)  of  Equation  (2)  when 

TL  .  2 

„/■  s  c  lOU 

f(u)  =  -  —  e 

where  f(u)  is  the  forcing  function  for  the  thrust  misalignment, 


We  get 


t  ,  2  X  /  TL  \ 

o  o 


2  2 
iau  -ioqu  .  - 

e  dudx  , 


TL  t  .  2  ^  2 

_  _  _c  j*  giqoQC  J^^ia(l-q)u 

^  t  t  dudx 

o  o 

Letting  A  =  oq  and  D  =  a(l-q),  then  Equation  (5)  reduces  to 
^  TL  t  ,,  2  X  -T,  2 

L(t)  -  -  -r  1-= 

t  t 

o  o 


Case  b) 


Dynamic  Unbalance 


We  let  equeil  the  ^  (t)  of  Equation  (2)  when 

f(u)  =  -  (l-2q)p^e^°^  (ia  -  (ot)^) 
where  f(u)  is  the  forcing  function  for  dynamic  unbalance. 
We  obtain 


^^(t)  =  J'  J  -(l-2q)S  (ia  -  (au)^)e'^°*^^  du  dx. 
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Again  letting  A  =  oq,  D  =  a(l-q),  and  rearranging  terms  we  get 

t 


=  -(l-2q)p^  J 


2  p  iAx^ 

-a  J  e 


2  X  2 

lAx  c  ^ 

e  J  ®  du  dx 

t 

0 


r,  2  IDu  ^  , 

J  u  e  du  dx 

t 


)  • 


(5) 


A  look  at  Equation  (4)  and  Equation  (5)  indicates  that  we  are 
faced  with  the  evaluation  of  the  following  integrals: 


t  .  2  X 

P^(t)  =  J-  .1*-  ; 

t  t 

o  o 


e  du  dx 


,  (6) 


and 


F2(^')  = 


t 

1 


iAx 


X 

1-^ 


■2  iDu^  -  , 

e  du  dx 


(7) 


Then  we  have 


TL 


and 


We 


0l(t)  =  -  ^  Fj^(t)  , 

yJ^Ct)  =  -  (l-2q)p^(iaF^(t)-a^F2(t)  )  . 

proceed  to  evaluate  F^(t)  and  F^Ct).  Define 


A 

■S' 


Fj(x)  =  \  e^^  du 


(8) 

(9) 

(10) 


Letting  x  =  ms  where  m  = 

x/m 


we  obtain 


F,  (.— )  =  m 

5  m 


1 


irtG^/2  , 
e  '  ds. 


'’"o^m 
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We  can  rewrite  Fj^(t)  as 


F^(t)  =  +  idj  +  3j^)) 

(14) 

where 

j  ..  pcos  Ax^(C(x/m)  -  C(t  /in))dx  , 

(15) 

^  t 

o 

t  o 

=  r  sin  Ax  (S(x/m)  -  S(t  /m))dx  , 

(16) 

^  t  ° 

o 

t 

=  Xcos  Ax^(S(x/m)  -  S(t^/m))dx  , 
t 

(17) 

0 

an4  t 

I,  =  X  sin  Ax^(C(x/ni)  -  C(t  /in))dx 
^  t  ° 

(18) 

0 

To  evaluate  F2(t)  as  defined  in  (7)  introduce 

the  notation 

F.  (x)  =  /  u^e^""  du  . 

(19) 

Letting  r  =  u,  dv  =  ue  du,  and  integrating  by  parts,  we 
obtain 

_  2  iDt  ^  .  X  2 

■n  /  \  1  /  iDx  ,  o  \  ,  1  f  iDu  , 

t 

o 


i 

2D 


(t^e  -  xe  )  +  m(c(x/m)  -  C(t^/ra)(20) 


+  iS(x/m)  -  iS(t^/rn))  ^ 
We  can  use  Fj^(x)  in  evaluating  F2(t). 

Fp(t)=/e^  F.(x)dx 
^  t 

o 


(21) 
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15 


1 

2i(A+D) 


i(A+D): 


(25) 


2i(A+D) 

1 

2(A+D) 


^  cos(A+D)t^ 
(^sin(A+D)t^ 


cos(A+D)t^^  ) 

sin(A+D)t  ^  ) 

o  J 


The  third  integral  on  the  ri^t  side  of  Equation  (2l)  is  equal  to 
t 


m 


t  2  2  r 

J  (cos  Ax  +  i  sin  Ax  )  I  C(x/m)  -  C(t^/m)  +  i(S(x/in)  (24) 
^o 

-  S(tym))j  dx  =  m(l^  -  +  idj  +  Ij^))  . 


Let 


=J 


^  cos  Ax^dx  , 


(25) 


and 


l6  =  J  sin 


Ax  dx 


(26) 


Combining  information  from  Equations  (2l),  (22),  (25), 
(24),  (25)  and  (26),  we  find 


cos  Dt  ^  +  it  cos  Dt  ^  +  it  sin  Dt  ^ 

o5o  000  o5 


^2^^^  2D  [^^o 

-  t^  sin  Dt^^  (cos(A+D)t^  -  cos(A+D)t^^) 

■  (si^(A+D)t^  -  sin(A+D)t^^)  +  m(l^  -  +  i(l^  +1^)) 

^o  t^2^  o.T>i.2T  1 

=  -  —  cos  Dt  I-r  “  sin  Dt  -  r— 

2D  o  6  2D  o  5  4D(A+D) 


(cos(A+D)t^  -  cos(A+D)t^^) 


(27) 
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l2(t) 


-  (l-2q)  I^)  -  °  t^(co3  Dt/  Ig 

2 

+  sin  Ij)  -  (coB(A+D)t^  -  oos(A+D)t^^) 

<^-^V  (29) 

r  n —  (i-2q)a^t  2 

+  1  [-(i-2q)a  (I^  -  Ij)  +  - 25— ^(=0=  »to  S 

-  sin  Ig)-  (sin(A+D)t^  -  sin(A+D)t^^) 


15 


A  METHOD  FOR  DETERMINING  EFFECTIVE  VALUES  OF 
DYNAMIC  UNBALANCE  AND  THRUST  MISALIGNMENT 

We  could  idealize  the  physical  setup  by  imagining  the  projectile  with 
its  center  of  gravity  somehow  held  fixed  in  space,  but  with  the  projectile 
free  to  move  about  its  center  of  gravity.  If  we  assume  that  the  transverse 
angiilar  motion  is  caused  entirely  by  thrust  misalignment  and  dynamic  un¬ 
balance,  solutions  to  Equation  (2)  can  be  used  to  provide  an  indirect 
method  of  estimating  L^  and  Indeed,  for  positive  unit  values  of  L^ 

and  6  we  could  obtain  the  solutions 
c 

-= -  and  — -  for  discrete  values  of  t  between  the  ignition  and 

Li  p 

c  c 

the  end  of  burning  times.  The  solutions  of  Equation  (2)  would  be  deter¬ 
mined  either  by  solving  it  directly  by  numerical  methods  or  from  Equations 
(28)  and  (29). 

We  would  then  obtain  an  over -determined  system  of  equations. 

0R(t^)  =  x(t^)Re(L^)  +  y(t^)Re(p^)  , 

0R(t2)  =  x(t2)Re(L^)  +  y(t2)Re(p^)  , 

(50) 


0R(t^)  =  x(t^)Re(L^)  +  y(t^)Re(p^) 

where 

0l(t) 

x(t)  =  real  component  of  — = - 


02(t) 


y(t)  =  real  component  of 


> 


Re(L^)  =  real  component  of 


Re(p^)  =  real  component  of 

As  mentioned  above,  x(t^)  and  y(t^)  are  obtained  from  the  solutions 

of  Eq^uation  (2).  The  angles  0(t^)  would  be  measured  experimentally 

(perhaps  by  a  transducer).  The  values  of  Re(L  )  and  Re(p  )  would  be 

c  c 

ualnovn.  Since  we  could  expect  errors  in  the  0,  x,  and  y  values,  estimates 

of  Re(L  )  and  Re(p  )  would  be  obtained  by  applying  least  squares  techniques, 
c  c 

Similar  expressions  would  provide  estimates  for  the  imaginary  components 

of  L  and  6  . 
c  ■  c 

IV  d. 

W.  J.  SACCO 
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